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Abstract 

A new parameter space is used to classify circular orbits in the Schwarzschild 
metric. 

PACS numbers: 04.20. Jb, 02.90.+p 

In Newtonian mechanics, the answer to the question "Is a circular orbit 
possible for a particle in a gravitational field?" is a simple yes. In general 
relativity, the answer that is given in most texts on this subject is not based on 
the exact solution of the orbit in the Schwarzschild metric, but is based on an 
analysis that makes use of the effective potential. Our treatment of this question 
is based on the explicit analytic expressions for the trajectories of a particle with 
a non-zero mass in the Schwarzschild metric that include bound, unbound and 
terminating orbits that we presented recently [1,2]. An understanding of how 
the various circular orbits given by the effective potential analysis are related 
to the exact solutions is important for clarity and for giving a definitive answer 
to the question. 

As we pointed out in ref.l, the trajectories of a non-zero mass particle can 
be conveniently characterized and placed on a " map" in a dimensionless param- 
eter space (e,s), where we have called e the energy parameter and s the field 
parameter. The parameter space (e, s), where < e < oo, and < s < oo, can 
be divided into two regions: Region I characterized by < s < Si, where si(e) 
is well defined, allows terminating orbits as well as the elliptic- (0 < e < 1), 
parabolic- (e = 1), and hyperbolic-type (e > 1) orbits, while Region II char- 
acterized by s > si allows terminating orbits only. It will be seen that this 
parameter space needs to be extended in order to represent all possible circular 
orbits. 

The Schwarzschild geometry is the static spherically symmetric gravitational 
field in the empty space surrounding some massive spherical object (which we 
shall refer to as a star or a black hole even though this also includes the case of 
a planet such as Jupiter or Earth) of mass M. We will consider a test particle 
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(which can be a spacecraft, an asteroid, or a planet), with its position relative 
to the star or black hole described in the spherical coordinates r, 0, <p, moving 
in the equatorial plane 6 = n/2. Let E be the total energy of the particle in its 
orbit and m its rest mass, and define 



moc- 
Let 

h = r 2 <f>, (2) 

which can be identified as the angular momentum per unit rest mass of the 
particle, and where the derivative ■ represents d/dr and r is the proper time 
along the path. The Schwarzschild radius a is defined by 

2GM 

a=—, (3) 

where G is the universal gravitation constant, and c is the speed of light. 
Making the substitution 

u=\, (4) 
the 'combined' energy equation [3] 

.2 - - - 2 

together with cq.(2) can be expressed as 



h 2 ( a\ c 2 a 9/ 9 „. 
+ ^( 1 --)- — = c2 ^ 2 - 1 ) ( 5 ) 



du\ 2 9 c 2 (k 2 - 1) c 2 a o 

As shown in rcf.l, it is more convenient to work with the dimensionless 
parameter U defined by 

in terms of which eq. (6) can be expressed as 

2 



dU 



= AU* - g 2 U - g 3 , (8) 



where 



1 2 

.92 = Y2 -s, 

93 = —-—s 2 + -(l-e 2 )s\ (9) 
y 216 12 4 v ; ' ' 
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and where the dimensionless parameter 



h 2 c 2 ( K 2 -iy 1/2 

i + 



that ranges from to oo is called the energy eccentricity parameter or simply 
the energy parameter, and the dimensionless parameter 

that ranges from to oo is called the gravitational field parameter or simply 
the field parameter. It is seen that the energy parameter e actually depends on 
the energy, the gravitational field, and the angular momentum of the particle, 
and the field parameter s depends on the gravitational field and the angular 
momentum of the particle, but not on the energy. 

The parameter space in our work presented in refs.l and 2 was confined to 
e > and s > 0. General relativity in fact allows an extended parameter space 
that includes e 2 < 0. Let us first note what e 2 > and e 2 < mean physically. 
From eqs.(2), (5), and (10), we find that e 2 > or < respectively imply 

.2 / • GmV 2GM 2 

r + I rq> > or < — ^ — r( P ■ 

\ r 2 <j> J c 

The right-hand side of the above inequalities is in the Newtonian limit, 
and thus it is important to note that while the parameter space e 2 > has its 
Newtonian correspondence, the "extended" space characterized by e 2 < does 
not have a Newtonian correspondence. In the following discussion, we shall first 
confine our results to the parameter space e > 0, and then discuss the additional 
results from those that include the extended parameter space. 

The discriminant A of the cubic equation 

4C/ 3 - g 2 U - g 3 = (12) 

is defined by 

A = 27 5 2 - g\ (13) 

and divides the parameter space (e, s) for e > and s > into two regions: 
Region I corresponding to A < covers the values of the field parameter < 
s < si, and Region II corresponding to A > covers the values of s > si, where 
Si as a function of e is given by [1] 



c2 _ 1 - 9e 2 + y/(l - 9e 2 ) 2 + 27e 2 (l - e 2 ) 2 
27(1 -e 2 

for e ^ 1, and s\ — 1/16 for e = 1. We note that the parameter value s\ 
that characterizes the (upper) boundary of Region I decreases very little from 
si = (2/27) 1 / 2 = 0.272165.. for e = to Sl = 1/4 = 0.250000 for e = 1, but that 
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it decreases more quickly as e increases beyond 1 and s\ — >■ (V27e) _1 — »■ as 
e — > oo [4]. Region I has terminating and non-terminating orbits that include 
bound (for < e < 1) and unbound orbits (for e > 1), and Region II has 
terminating orbits only. 

In Region I (A < 0), the three roots of the cubic equation (13) are all real. 
We call the three roots e\, e 2 , e 3 and order them so that e\ > e 2 > e 3 . Define a 
dimensionless distance q measured in units of the Schwarzschild radius a by 

q = r ~, (15) 
which is related to U of eq.(7) by 

I = i+4tt (16) 

The equation for the non-terminating orbits in Region I, bound and un- 
bound, and for < e < oo, is given in terms of Jacobian elliptic functions [5] of 
modulus k by 

- = i+4e 3 + 4(e 2 -e 3 )sn 2 ( 7( /.,fc), (17) 
q o 

where e\ > e 2 > U > e 3 [6]. The constant 7 appearing in the argument and 
the modulus k of the Jacobian elliptic functions are given in terms of the three 
roots of the cubic equation (12) by 



ei - e 3 

and the roots e\, ei, e 3 are given by 



7 = (ei-ea) 1 / 2 , 

e = (18) 




where 

/27\ 1/2 

cos6 = g 3 ^j . (20) 

The modulus k of the elliptic functions has a range < k 2 < 1 as 9 varies 
from to 7r. The period of sn 2 (-f<f),k) is 2K(k), where K(k) is the complete 
elliptic integral of the first kind [5]. For k = 0, sn(x, 0) = sin a;, cn(x, 0) = cosx, 
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dn(x,0) = 1, and for k 2 = 1, sn(70, 1) = tanh(70), cn(70, 1) = dn^cj), 1) = 
sec/i(7^>). As fc 2 increases from to 1, i4T(fc) increases from 7r/2 to oo. 

The orbits given by eq.(17) for < e < 1, e = 1, and e > 1 will be referred to 
as the elliptic-, parabolic- and hyperbolic- type orbits respectively, even though 
their shapes can differ greatly from the conic sections that arise in Newtonian 
mechanics, as shown in ref.l. They correspond to the cases for which < k 2 < 1, 
k 2 = 1 , and k 2 > 1 respectively, where k is defined by eq.(l). 

For the elliptic-type (0 < e < 1) orbit, the distance r of the planet from the 
center of the star or black hole assumes the same value when its polar angle <j> 
increases from (f> to <j) + 2K/j. Comparing this with the increase of 4> from </> to 
(f> + 2tt in one revolution for the planet, the precession angle is the difference of 
the two and is given by 

A* -2*. (21) 

7 

For k 2 close to the value 1, the particle can make many revolutions around 
the star or black hole before assuming a distance equal to its initial distance. 

The maximum distance r max (for < e < 1) of the particle from the star 
or black hole and the minimum distance r m ; n (for < e < oo) of the particle 
from the star or black hole, or their corresponding dimcnsionless forms q max 
(= r mM /a) and q nlin (= r min /a), are obtained from eq.(17) when 70 = and 
when 70 = K(k) respectively, and they are given by 

— = \+^z, (22) 



and 
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q +4e 2 . (23) 

In contrast to the energy eccentricity e defined by eq.(10), the true eccentric- 
ity e of the elliptic- and parabolic-type type (0 < e < 1) orbit given by eq.(17) 
is defined by 

''max ^min Qmax Qmin &2 63 . . 

£= T = 7 = TTTZ — T — ' ( 24 ) 

'"max T" ^min <?max + <?min 1/0 + e 2 + e 3 

and the true eccentricity e coincides with the eccentricity e only when s = for 
0<e<l(e = l when e = 1 for all values of s). The relationship between e 
and s was exhibited in ref.l. 

We now describe in some detail two special cases (i) k 2 ~ and (ii) k 2 = 1 
which relate to the possibility of circular orbits for the parameter space being 
considered (e > 0). 

(i) The special case of k 2 ~ 0. 

First we note that if k 2 = exactly, then it follows from eqs.(18) and (19) 
that this requires 6 — 0, and thus e\ = 2(g 2 /12) 1 / 2 , e 2 = e 3 = -(52/I2) 1 / 2 , and 
from eq.(20) that s = 0, from eq.(13) that A = 0, and from eq.(24) that e = 0. 
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Thus wc find e\ = 1/6, = = —1/12, and 7 = 1/2. Because s = implies 
zero gravitational field, it is important to recognize that the classical Newtonian 
case corresponds to the case of s ~ and k 2 ~ 0, and not s — and k 2 = 
exactly. 

Consider the Newtonian limit s ~ and k 2 ~ 0. The constant c 2 (k 2 — 1) 
which is < for a bound orbit and > for an unbound orbit, can be identified 
with 2Eo/m in the Newtonian limit, where Eq is the sum of the kinetic and 
potential energies and is < for a bound orbit, and is > for an unbound 
orbit, and m is the mass of the particle (which approaches mo). In this limit 
we therefore obtain 



for e > 0, while the case for e = will be analyzed more carefully below. 
We note that, just like the Newtonian orbit equation above for which the cases 
< e < 1 and e > 1 give rise to bound (elliptic) and unbound (parabolic and 
hyperbolic) orbits, the single analytic expression represented by eq.(17) with the 
same criteria for e gives rise to precessing bound (elliptic-type) and unbound 
(parabolic- and hyperbolic-type) orbits. We should also note that in general 
relativity a second parameter s with a range < s < s±, where s\ is given by 
eq.(14), is needed, for these three types of orbits to be realized. For the case 
s > si, all orbits are terminating. 

We now consider the Newtonian limit characterized by very small values of 
s and k 2 for the case e = more precisely. Setting e = in g 3 in eq.(9) and 
using eq.(20), the expansion of cos6> is found to be 




1 GM 



(1 — ecos(/>) 



costf = 1 - 2 2 • 3 3 s 6 - 2 • 3 5 • 5s 
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(25) 



from which wc find 




(26) 



From eq.(19), we then calculate the expansions [7] 



d = \ -s 2 -3s 4 -20s 6 + ... 
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s 5 + 10s 6 + .... 



s 5 + 10s 6 + ... 



(27) 
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From eq.(18), we obtain 



1 3 2 V2 , 27 4 33V2 - 

k 2 - 4V2.S 3 (l + ^s 2 + , (28) 
and from [5] K{k) = (tt/2)(1 + k 2 /A + ...), we find 

K (k) = |(l + V2s 3 + ...) . (29) 
Using the above equations and eq.(19), the orbit equation (17) becomes 

i ~ 2s 2 |l - V2.s [l - 2.sn 2 (7</-, fc)] | , 
or, after substituting eqs.(3), (11) and (15), 

-^^{l-V2 S [l-2,sn 2 (70,fc)]}, (30) 

where 7 and A; 2 are given by eq.(28). The elliptic function sn(^<t>, k) can be 
expanded in power series in s for e = as 

(a n • ^ 3< ^ 2 . 
sn(j(p 7 k) — sin — — s cos — + ... 

To order s, we find that the Newtonian limit of the orbit equation for e = 
can be expressed as 

I = ^(l-^ SC os[(l-^]), (31) 

where 6 ~ 3s 2 . Thus, according to general relativity, the Newtonian limit 
(for s and k 2 very small) of the orbit of a particle for e = is not a true circular 
orbit but is an orbit with an eccentricity equal to y/2s = V2GM/(hc), and is 
thus an elliptical orbit. 

Using eqs.(24) and (27), we find the true eccentricity is given by 

e = Vis + ... 

in agreement with that given in eq.(31). Only if s = does the orbit rep- 
resented by eq.(31) become a true circular orbit, but then it is the case of zero 
gravitational field and the radius h 2 /(GM) of the orbit is infinite (h 2 /(GM) = 
(hc/GM)(h/c) — (oo)(/j/c) because s = GM/(hc) = 0), and thus the particle 
is going in a straight line. 

The precessing elliptical orbit for small s is represented by eq.(31) or eq.(30). 
The angle of precession, from eq.(21), is given by 

. t „ 2 6tt(GM) 2 

A^6ns 2 = ^-^. (32) 
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We recall from previous work [1, 2, 7] that for e > 0, 7 = [1 — (3 — e)s 2 + ...], 
k 2 = 4es 2 + the approximate precessing elliptical orbit equation is 1/r = 
{GM/h 2 ){\ - ecos[(l - 5)(f>]} or 1/r = (GM/h 2 ){l - e[l - 2sn 2 ( 7 </>, jfe)]}. Note 
that to the order of s 2 , we get the same precession angle A(j> given by eq.(32) 
for all elliptical orbits independent of e. 

We now consider the second special case for k 2 . 

(ii) The special case of k 2 = 1. 

It follows from eqs.(18)-(20) and (13) that cos# = —1, e\ = e 2 = ^Jg^/VZ, 

e 3 = — \/g2/3, 7 = (3g2/4) 1/ ' 4 , A = 0. The range of s values is given by 
eq.(14) and gives the boundary between Regions I and II. On this boundary 
dchncd by k 2 — 1 we have 32 > 0, 33 < 0, and the relation between the two is 
^93 = -1/52/3. 

For k 2 = 1, eq.(17) becomes 



q 3 ' "V 12 - Vl2 sec M7<W (33) 

The orbit will be referred to as an asymptotic orbit. For < e < 1, the 
particle starts from an initial position <7 max at <j> = and ends up at </> = 00 
circling the star or black hole with a radius that asymptotically approaches q m in, 
where q mm is given by 

where g 2 is evaluated using s — s x given by eq.(14). For e > 1, the particle 
starts from an initial position at infinity at an angle ^1 to the horizontal axis 
given by 

2, - v 2 1 :5 



tanh (7 * l) = 3 18V52 

and ends up at <j> = 00 circling the star or black hole with a radius that 
asymptotically approaches q m i n - 

A possible case of circular orbits is the case of those limiting circles that the 
asymptotic orbits given by eq.(33) become. The range of the radius q c given by 
9min of eq.(34), using the values of si given after eq.(14), is between 2.25 for 
e = and 1.5 for e = 00, i.e. the radii r c of the limiting circular orbits are in 
the range 

3GM 4.5GM 

<r c < — . 35 

and they occur for the values of s = S\ given by eq.(14). However, it is clear 
that these circular orbits are unstable because orbits exist which are spirally 
asymptotic to them [8]. 

Another possible example of a circular orbit comes up as a special case of 
the terminating orbits in Region I. The terminating orbit in Region I is given 

by [1] 
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1 = 1 , fi ~e 2 sn 2 (7^, fc) 

q 3 cn 2 ( 7 (/.,fc) ' 1 J 

where [/ > ei > e2 > e3 and 7, fc, ei, e2, e 3 are given by eqs.(18) and (19). 
For < e < 1, the particle, starting from the polar angle <j) = in a direction 
perpendicular to the line joining it to the star or black hole, at a distance q\ 
from the star or black hole, plunges into the center of the black hole when its 
polar angle reaches 4> 1 , where qi and (f> 1 are given by 

- = ^+4e l5 (37) 
qi o 



and 



K(k) 



1 

For e > 1, the particle starts from infinity at a polar angle (f> = ^1 > to 
the horizontal axis given by 

sn 2 7 *! = -77 r , 

4(e 2 - e 3 ) 

and plunges into the center of the black hole when its polar angle reaches 

K 

The terminating orbit given by eq.(36) becomes a circular orbit when k 2 = 1, 
because for k 2 = 1, the two roots ei and e 2 become equal (see the description 
before eq.(33)) and 1 — tanh 2 (7</>) = sech 2 ("/4>). We find that the radius q c of 
the circular orbit is given by 

H+M+VI- <38 > 

which is the same as the radius of the limiting circular orbit given by eq.(34). 
Thus q c ranges from 2.25 for e = to 2 for e = 1 and to 1.5 for e = 00, and 
the radius r c of the circular orbit is in the same range as that given by eq.(35). 
However, it is obvious that these circular orbits are unstable because a small 
perturbation would make them into terminal orbits. 

We now consider the relationship of these results from the exact solutions 
given by eqs.(17) and (36), in particular eqs.(22), (23), (34) and (38), to the 
solutions for the circular orbits obtained from the analysis using the effective 
potential [9], or from the simple consideration of r — const.. Differentiating 
eq.(6) with respect to 0, we obtain 

d 2 u c 2 a 3a 9 

w +u =^ + ^ u - (39) 

Setting d 2 u/d(j) 2 = gives a quadratic equation in u and the two solutions 
for u can be expressed in q given by eq.(15) and in g 2 given by eq.(9) as 
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1 1 ±4 /»,N'/ 2 



(ft) («> 
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The effective potential analysis [9] chose what seemed to be a reasonable 
range 

< s 2 < y 2 (41) 

and thus 

^ > .92 > (42) 

for the solutions given by eq.(40) to remain real. Applying the inequalities 
given by eq.(42) to the solution given by eq.(40) with the lower sign (corre- 
sponding to the minimum of the effective potential) gives a range 3 < q < oo or 
a range for the radius r of a possible stable circular orbit as 

—5- <r<oo, 43 
cr 

and applying the same inequalities given by eq.(42) to the solution given 
by eq.(40) with the upper sign (corresponding to the maximum of the effective 
potential) gives a range 1.5 < q < 3 or a range for the radius r of a possible 
unstable circular orbit as 

3GM 6GM 

— (44) 

We now want to see what our analytic solutions give for the parameter space 
we are using. 

First we realize that eq.(17) represents a circular orbit only when &2 = 
e^. Because is negative for the entire range < k 2 < 1 corresponding to 
< 8 < 7r (see eq.(19)), q must be greater than 3 and hence there is a lower 
bound to the radius r > 6GM/c 2 which is in agreement with eq.(43). However, 
equating the expressions for ei and given in eq. (19) requires that 8 = if 
we confine ourselves to e > so that we have €2 = £3 = -(32/12) 1 / 2 . Under 
these conditions, one has, from eqs.(22) and (23), 

ymin ymax ° <-> 

which can be identified with the solution (40) with the lower sign. However, 
setting e2 = e% in eq.(18) shows that the solution given by eq.(45) requires 
k 2 = which, together with 9 = also requires s = as the only possible 
value for s (which gives r = 00) (see the discussion of the special case of k 2 ~ 
after eq.(24) and the discussion after eq.(31)). Thus our exact solution excludes 
any finite value of r given in eq.(43) as a possible radius for a circular orbit 
for the parameter space we are using. But s = corresponds to the case of 
zero gravitational field which makes this solution useless. Our solution given 
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by eq.(30) gives the best possible nearly circular orbit; it states that a stable 
most "circular" orbit with e = and with s ~ and k 2 ~ must have a small 
but nonzero true eccentricity e equal to y/2s and is thus elliptical and precesses 
with a precession angle A(j> ~ 67rs 2 . 

We now consider the solution given by the upper sign in eq.(40). It can be 
identified with our analytic solutions (17) and (36) when k 2 — 1 (see eqs.(34) 
and (38)). Note that these are two separate cases: Eq.(34) gives the limiting 
radius of an asymptotic orbit from solution (17), and eq.(38) gives the radius of 
a circular orbit that a terminating orbit given by solution (36) becomes under a 
special case. As previously discussed, our exact solutions clearly show that both 
these cases correspond to unstable circular orbits. Our exact analytic solutions 
that apply for A < show the range given by eq.(35) to be the possible range 
of radii for an unstable circular orbit for the parameter space we are using. 

Thus comparing our analytic solutions while confining ourselves to the pa- 
rameter space e > with those obtained from the effective potential theory, 
we find a stable circular orbit only for r — oo compared with eq.(43), and an 
unstable circular orbit in the range given by eq.(35) compared to that given by 
eq.(44). 

We now consider the extended space characterized by e 2 < 0. Our analytic 
solution (17) still applies in this region given by A < 0. The upper boundary of 
this region is still given by eq.(14) from the coordinate point e 2 = 0, s 2 = 2/27 
to the point e 2 = —1/3, s 2 = 1/12, which we shall call the vertex point, at 
which <?2 = .93 = 0, and the three roots of the cubic equation (12) are equal to 
zero. The upper boundary characterized by k 2 = 1 in the extended parameter 
space is still defined by si (eq.(14) with 2/27 < s\ < 1/12) and eqs.(34) and 
(38) still hold. Thus in the extended space we have unstable circular orbits with 
radii in the range given by 

4.5GM 6GM 

— < r ^ 46 

c C 

which, together with the possible unstable circular orbits with radii given 
by eq.(35) in the parameter space e > 0, make up the total interval given by 
eq.(44) which was derived by the effective potential theory. A similar analysis 
holds for the range of radii of the limiting circles of the asymptotic orbits. 

The left boundary of Region I is given by s = s[ , where s[ is given by 

a _ l-9e 2 -y/(l-9 e 2 ) 2 + 27 e 2 (l- e T 
Sl " 27(1 -e 2 ) 2 ' [ ' 

and it extends from the vertex at (e 2 ,s 2 ) = (—1/3,1/12) to the origin at 
(e 2 ,s 2 ) = (0,0). It can be verified that along this left boundary of Region I, 
cos 9 given by eq.(20) is equal to 1, and we have ei = e$, k 2 = 0, and e = from 
eq.(24). This means that we have possible stable circular orbits along the left 
boundary curve with radii in the range given by eq.(43). It is to be noted that all 
possible stable circular orbits with finite radius are in the extended parameter 
space e 2 < 0, and the physical implication of this extended parameter space 
should be noted. 
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In summary, the analysis of our exact analytic solutions for a non-zero mass 
particle in a gravitational field in the Schwarzschild geometry shows the follow- 
ing results: 

(1) A stable circular orbit is not possible in the parameter space e > 0. The 
most " circular" orbit of a particle in a gravitational field must have a small but 
non-zero eccentricity and is thus an elliptical orbit and it precesses. For e = 
and small values of s and k 2 , the orbit equation is given by eq.(30). However, 
in the parameter space e 2 < 0, a stable circular orbit with a radius given by 
eq.(45), where s = s^, and with a range given by 6GM/c 2 < r < oo, is possible. 

(2) An unstable circular orbit with a radius given by eq.(38), where s = 
si, and with a range given by 3GM/c 2 < r < 4.5GM/c 2 , is possible in the 
parameter space e > 0, and with a range given by 4.5GM/c 2 < r < QGM/c 2 in 
the parameter space e 2 < 0, . 

We will present a more complete analysis in future work. 
References 

^Electronic address: fhioe@sjfc.edu 

[1] F.T. Hioe and D. Kuebel, Phys. Rev. D 81 , 084017 (2010). 

[2] F.T. Hioe and D. Kuebel. larXiv:1008.1964l vl (2010). 

[3] M.P. Hobson, G. Efstathiou and A.N. Lasenby: General Relativity, Cam- 
bridge University Press, 2006, Chapters 9 and 10. 

[4] The misprint s\ — > (v^e) -1 on the 5th line of Section V in ref.l should 
be corrected to s\ — > (\/27e) _1 . 

[5] P.F. Byrd and M.D. Friedman: Handbook of Elliptic Integrals for Engi- 
neers and Scientists, 2nd Edition, Springer- Verlag, New York, 1971. 

[6] Equations similar to eq.(17) in somewhat different forms were given by 
A.R. Forsyth, Proc. Roy. Soc. Lond. A97, 145 (1920), C. Darwin, Proc. Roy. 
Soc. Lond. A249, 180 (1958), ibid. A263, 39 (1961) and S. Chandrasekhar: 
The Mathematical Theory of Black Holes, Clarendon Press, Oxford 1992, p. 
106-109. 

[7] The coefficients beyond s 2 in e2 and in eq.(20) in F.T. Hioe, Phys. 
Lett. A 373, 1506 (2009) are incorrect. However those errors did not affect any 
other results presented in that paper. 

[8] These orbits were given in E.T. Whittaker: Analytical Dynamics, 4th 
Edition, Dover Publications 1937, p. 393. 

[9] See e. g. ref. 3 Sections 9.8 and 9.9, also J.V. Narlikar: An Introduction 
to Relativity, Cambridge University Press, 2010, Section 9.3. 



12 



